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ABSTRACT 

High-quality velocity maps of galaxies frequently exhibit signatures of non-circular 
streaming motions. We here apply the software tool, velfit recently proposed by 
Spekkens & Sellwood, to five representative galaxies from the THINGS sample. We 
describe the strengths and weaknesses of the tool, and show that it is both more 
powerful and yields results that are more easily interpreted than the commonly used 
procedure. We demonstrate that it can estimate the magnitudes of forced non-circular 
motions over a broad range of bar strengths from a strongly barred galaxy, through 
cases of mild bar-like distortions to placing bounds on the shapes of halos in galaxies 
having extended rotation curves. We identify mild oval distortions in the inner parts 
of two dwarf galaxies, NGC 2976 and NGC 7793, and show that the true strength of 
the non-axisymmetric gas flow in the strongly barred galaxy NGC 2903 is revealed 
more clearly in our fit to an optical Ha map than to the neutral hydrogen data. The 
method can also yield a direct estimate of the ellipticity of a slowly-rotating potential 
distortion in the fiat part of a rotation curve, and we use our results to place tight 
bounds on the possible ellipticity of the outer halos of NGC 3198 and NGC 2403. 

Key words: galaxies: haloes - galaxies: kinematics and dynamics - galaxies: spiral 
- galaxies: structure - methods: data analysis 



1 INTRODUCTION 



The centrifugal balance of gas on near circular orbits in 
a galaxy yields a direct estimate of the central attraction, 
which is the first step towards a model for the mass distri- 
bution. The estimatio n of galaxy rotation c urves therefore 
has a long history (see lSofue fc Rubinll200ll . for a review). 

Well-sampled 2D velocity maps provide sufficient infor- 
mation to identify the rotation centre and to test whether 
the line-of-sight velocity field is, or is not, consistent with 
a circular flow patttern in an inclined plane about a com- 
mon rotation ce ntre. The widely-used software utility rotcur 
|Begemanlli987h divides the velocity map into a number of 
elliptical elements that are assumed to be projected circles 
around which the gas moves on circular orbits. It yields an 
estimate of the circular speed in each annulus, and has op- 
tions to fit for, or hold fixed, the rotation centre, systemic 
velocity, position angle and inclination in each annulus. It is 
uniquely powerful in its ability to extract information when 
the gas layer is warped. 

A number of possible systematic errors in the 
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fitted velocity have been d iscussed. Beam smearing 
(van den Bosch & Swatersll200ll ) is obviously reduced by im- 
proved spatial resolution of the observations, e.g., by using 
optical data when available. Various forms of turbulence - 
also described as pressu re support or as an asymmetric drift 
ijValenzuela et all\200i[ ) - can be recog nized and corrected 
for in high quality data l|Oh et al\ \2008). 

The present paper addresses the systematic error caused 
by non-circular flow patterns. Velocity "wiggles" or larger- 
scale distortions of the isovel ocity contours have long been 
recognized (e.g. 



Bosnia! Il978l) in wel l-sam pled 2D velocity 



maps. Both lHavashi fc Navarro! l|2006l ) and lValenzuela et al\ 
(2007) show qualitatively that non-circular streaming in 
non-axisymmetric potentials can lea d to a mis-est i matio n 
of the central attraction. Note that Weiner et ai. J200lf) . 



iKranz et all (|2003l ). iPerez et al] j|2004l ) and IWeinerl ([2004 1 



had previously published a number of quantitative models 
of non-circular flows that took this effect into account; these 
papers used the additional information contained in the non- 
circular motions of strongly barred or spiral galaxies to sep- 
arate the central attractions of the disc and halo. 

Not only do non-circular motions arise from bars and 
other visible distortions, but they may also be caused by ex- 
pected asphericities in dark matter halos (e.g. Ijing fc Sutol 
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120021 ; lAllgood et aZj|2006l ; lHavashi et al.ll2007l ). Because ha- 
los are largely pressure-supported, the principal axes of such 
distortions can slew only slowly with radius, and should ap- 
pear approximately straight over a moderate radial range. 
Even though halo shapes are expected to be made rounder 
by d isk formation (e.g. iDubinskH 1 19941 ; iDebattista et all 
2008), gas in the outer discs is predicted to flow in an el- 
liptical streaming p attern in response t o forcin g by a non- 
axisymmetric halo. Ha vashi fc Navarrol l|2006l ) argue that 
halo-driven non-axisymmetric streaming motions in the in- 
ner disc may mask a central cusp in the halo density. 

It is therefore desirable to be able to identify non- 
circular streaming motions in high-quality 2D velocity 
maps and to correct for their influence when deriving an 
estimate of the average c entral attract i on. T he power- 
ful approach developed by IWeiner et all {2001) (see also 
IZanmar Sanchez et q/.|[2008h requires multiple datasets and 
is time-consuming to implement. Thus there is a clear need 
for versatile tool for routine use. 

The procedure propo sed bv lSchoenmakers et all (|l997l ) 
and ISchoenmakera (Il999l ) estimates the magnitude of non- 
circular motions, but does not correct for them. It is 
embodied in the tool reswri, which is an extension of 
rotcur, that has been used quite extensively. More recently 
ISpekkens fc Sellwood! l|2007l . hereafter SS07) proposed an al- 
ternative tool, velfit, that does provide an improved estimate 
of the mean orbital speed. Here we demonstrate the superior 
performance of velfit by direct comparison of the results from 
several galaxies. We have chosen to ma ke this comparison 
using a few galaxies from the THINGS l|Walter et ai.ll2008l ) 
sample, for which non-circ ular motions have r e cently been 
estimated using reswri by iTrachternach et all (|2008l . here- 
after TBWBK). In addition, we show that velfit can also be 
applied t o optical data from a st rongly barred galaxy in the 
BHaBar (|Hernandez et <aZ.ll2005h galaxy sample. 



2 ESTIMATING NON-CIRCULAR MOTIONS 

In this section, we describe and compare the two separate 
software tools that are available for estimating the magni- 
tudes of non-circular flows. 

The procedure proposed bv lSchoenmakers et all (|l997f ). 
embodied in the tool reswri, is an extension of rotcur to 
include an harmonic analysis of the line-of-sight velocities 
around each ring. These authors use epicycle theory to in- 
terpret the fitted non-axisymmetric coefficients, and relate 
small values to the magnitude of the potential distortion. 

It should be noted that non-circular motions are readily 
confused with the kinematic signature of a warp, since both 
cause variations in the ellipticity and position angle of the 
flow pattern^ Thus, if the projection geometry is allowed 
to vary from ring to ring, then a large part of the actual 
non-circular motion may be masked by radial variations in 
the p osition angle (PA ) and inclination (i). Thus most users 
(e.g. iFathi et all 120051 ; Ivan Evmeren et aZ.| [2009) advocate 



1 In fact, the kinematic signatures of an oval potential and a 
warp are degenerate only when the principal axis of the poten- 
tial c oincides with either the major or minor axis of the projec- 
tion {Pence fc Blackman 1984; Fran x. van Gorkom fc de Zeeuwl 
Il994h . 



constraining i & PA to have the same values at all radii. 
On the other hand, TBWBK justified allowing individual 
tilts for rings in all the THINGS galaxies from the fact that 
the magnitudes of non-circular streaming motions in two 
galaxies (NGC 3198 & DDO 154) were little changed when 
i & PA were allowed to vary compared with when they were 
constrained to be constant with radius. Here we show that 
allowing rings to tilt independently led them to miss non- 
axisymmetric distortions in at least two galaxies: NGC 2976 
& NGC 7793. 

The tool reswri has a number of disadvantages. 

(i) A cos(m0) component in the projected velocities can 
arise from both cos[(m ± 1) 6] distortions to the poten- 
tial (ISchoenmakers et aZ.lll997l . see also SS07 and references 
therein), complicating the interpretation of the results. 

(ii) It assumes the perturbed velocities are small, and 
is therefore inadequate to characterize strongly non- 
axisymmetric flows. 

(iii) It makes no correction for the presence of non- 
circular motions. As a result, the estimated circular speed, 
V C (R), is biased toward the value of the azimuthal speed on 
the major axis of the projection, as explained below. 

(iv) It treats each ring independently, implying that a 
mild distortion that is coherent over a significant radial 
range is more easily masked by noise. 

(v) An estimate of the strength of the mild potential dis- 
tortion responsible for the detected non-axisymmetric flow 
requires an undesirable difference between two of the fit- 
ted coefficients, a consequence of point (i) above, and also 
includes the sine of an angle whose value cannot be deter- 
mined from this approach. 

SS07 therefore proposed an alternative tool, velfit illus- 
trated in Fig. [lPI It differs at root by postulating a spe- 
cific model of the flow that includes a possible non-circular 
streaming pattern about a fixed direction in the disc plane, 
i.e., a distortion that has no spirality. 

In particular, SS07 assume a flat disc plane and a dis- 
tortion having a fixed orientation at all radii, although its 
amplitude may vary with R. The tool velfit then fits the 
projected model velocity at a general point (eq. 5 of SS07): 



V S ys + sin i \Vt cos 9 

— V m ,t cos m(6 — 4> b ) cos 6 

— V m ,r sin m(8 — <j> b ) sin 6 ] , 



(1) 



to derive Vt(R), V m ,t(R), V m ,r(R), the angle 4> b , the systemic 
velocity V sys , and i & PA. (The various quantities are de- 
fined in the caption to Fig. [1] which illustrates an elliptical 
streaming pattern for which m = 2.) 

The code fits a model by minimizing the usual function 




(2) 



Here V b s (x,y) and a n are respectively the value of the 
observed velocity and its uncertainty for the n-th pixel at 
the position (x, y) on the sky. The index k ranges over K, 
which is the total all three sets of velocities, 14 , that define 

2 The velfit software is publicly available from 
http: / /www. physics . rutgers . edu/~spekkens/velf it/ 
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Figu re 1. A schema t ic illu stration of the approach taken by 
Spckk (.■lis fc Sellwoodl (20071) . The magenta lines in the upper 
panel indicate non-circular streamlines in the disc plane, which 
is shown face-on. The galaxy is observed in projection, with the 
intersection of the sky-plane and the disc plane making an angle 
(j>l to the major axis of the elliptical streamlines. Our model of 
the flow consists of a set of values for Vt, V2,t and V2,r around 
rings (dotted) at fixed radii. We project the model and fit to the 
observed data, using linear interpolation to predict data values 
between rings. 

The lower panels illustrate the bisymmetric variations of the ra- 
dial and azimuthal model velocities (arbitrary scale) in the disc 
plane around the cyan circle in the upper panel. The radial ve- 
locity (blue) varies with angle 6>j, = 6 — 4>b t° the major axis of 
the elliptical streaming pattern as — V2, r sin 2d^. The azimuthal 
velocity (green) varies as — Vn cos2$(, about the mean (Vt shown 
in red). This physically-motivated phase difference between the 
two non-axisymmetric components is embodied in the fit. 



the model velocity (eq. [TJ around each of the ellipses. The 
weights Wk, n , which include the trigonometric factors, also 
define the interpolation scheme that yields a model predic- 
tion at the projected position (x,y). The code allows the 
radial extent of the distortion to be restricted, if desired, 
while a simple axisymmetric flow is fitted over the remain- 
der of the data. SS07 describe the algorithm in detail. 

We here use the term "bar model" to denote any 
straight bisymmetric distortion no matter what its origin 
or amplitude. Fig. [1] should not, however, be interpreted to 
imply that the code can fit only bars, or oval features, al- 
though that is its most useful application. The same code 



can be used to fit models having higher or lower rotational 
symmetry, m, or axisymmetric (m = 0) radial flows. 

This tool avoi ds the above-listed d i sadva ntages of the 
method devised bv lSchoenmakers et al\ (|l997l ). as follows. 

(i) Fitting a specific distorted flow model to the projected 
data avoids the complications caused by the coupling of dif- 
ferent angular periodicities. 

(ii) It can fit for arbitrarily large distortions because it 
does not require V m ,t(R) & V m ,r(R) to be related by the 
epicycle approximation. 

(iii) It yields Vt, which is an improved estimate of the 
average orbital speed at each radius, as discussed below. 

(iv) It uses all the data in a single fit, making it eas- 
ier to identify coherent mild distortions in noisy data and 
to go some way towards "averaging over" small-scale spiral 
streaming. 

(v) The magnitude of a mild potential distortion is much 
more directly related to the fitted velocity coefficients, as we 
show in section [3] 

An elliptical flow pattern is clearly an idealization. It is a 
first-order improvement that captures the essential large- 
scale features of a bar-like flow, and will therefore be a better 
fit than a flat, axisymmetric model. It does not attempt to 
include other features, such as shocks in the bar, spiral arms, 
turbulence, warps, etc., that are also generally present in the 
data. 

The assumption of a flat disc plane is equivalent to re- 
quiring constant i & PA in reswri, which is the usual prac- 
tice, except that velfit has the further advantage that it de- 
termines the optimal values as part of the fit. The inner 
parts of spiral galaxy discs are believed to be flat: warps 
are gener ally observed to start near the edge of the optical 
disc (e.g. lBriggslll990l ) and theoretically we expect the mas- 
sive inner disk to be co herent enough to resist bending (e.g. 
IShen fc Sellwoodll2006h . 

With the projection angles held at constant values for 
all radii, it may seem from Fig. [T]that fitting a simple circu- 
lar flow model, e.g. with rotcur, would yield the same Vt as a 
bisymmetric model from velfit. However, the estimated cir- 
cular speeds for an axisymmetric model are biased towards 
values on the major axis where the line-of-sight component 
of the orbital motion is greatest. Since gas moving on el- 
liptical orbits has its smallest orbital speed at apocentre(f| 
the fitted Vt is biased low when the streaming pattern is 
closely aligned with the major axis. Conversely, we should 
expect an axisymmetric fit to be biased high by the higher- 
than-average speed of the gas near pericentre when the bar 
is oriented close to the projected minor axis, and to yield a 
fair estimate of Vt when the bar is at intermediate angles. 
Recall that Vt is defined to be the average tangential ve- 
locity around a circle that may cross multiple streamlines, 
as shown in the strongly non-axismmetric flow sketched in 
the upper panel of Fig. [T] Thus even though Vt affords a 
fairer estimate of the central attraction than that from an 



3 Even in strong bars, where angular momentum is not approx- 
imately conserved, gas at apocentre of the stream lines must be 
moving slower than the average since it plunges inwards after that 
point, irrespective of the angular rotational speed of the potential 
that drives the non-circular flow. 
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axisymmetric fit, it is not the circular speed in the equiv- 
alent axisymmetric potential, unless the distortion is very 
mild. 

Because reswri fits each ring independently, it is not 
easy to apply a smoothing constraint to the fit. Spiral arms, 
turbulence, etc., produce localized distortions to the flow 
that can lead to "wiggles" in the fitted velocities, as well as 
small-scale variations in V S ys, i, & PA. The user of reswri can, 
and probably should, hold these global parameters fixed, 
but the tool does not have the option to smooth the fitted 
velocities. While excessive smoothing is clearly dangerous, 
e.g. it could reduce the slope of the inner rotation curve, a 
small degree of smoothing applied to constrain the fit can 
be beneficial. In Appendix A we describe how an optional 
penalty can be applied within velfit to smooth the radial 
variations of the fitted functions Vt(R), Vm,t{R), & V m ,r(R)', 
the magnitude of the penalty can be set independently for 
the axisymmetric and non-axisymmetric terms. 

The tool velfit assumes a straight position angle for 
the non-axisymmetric distortion, which has both advantages 
and disadvantages. Since it cannot follow the radial winding 
of a spiral pattern, it "averages over" the spirals (advan- 
tage (iv) mentioned above) and fits only straight distortions 
with m-fold rotational symmetry. Its primary purpose there- 
fore is to identify bar-like or oval distortions in the disk or 
halo. However, a fixed position angle precludes mapping the 
distorted flow of a spiral, which could be an interesting ca- 
pability. It may be possible to adapt velfit to include extra 
parameters to fit a spiral of a certain shape, which we leave 
for possible future development. 

A weakness of velfit is that it tends to return absurd 
velocities when the bar angle, <j>t, is near zero or 90° be- 
cause a degeneracy arises between the velocity components 
at these special orientations. To see this, consider equation 
{TJ for the predicted projected velocity when cj>i, = 0: the 
products cos 29 cos 6 and sin 26 cos 9, can be separated into 
a part that varies as cos 9 and another that varies as cos 39 or 
sin 39. Therefore Vt is partly degenerate with both Vij, and 
V2, r , and a similar partial degeneracy arises when 4>b = 90°. 
In principle, the 39 variation of the model breaks the de- 
generacy, but these more rapid angular variations are more 
susceptible to noise, and the "best fit" values of the three 
velocity components can be absurd. We show an example in 
section |4~1I where smoothing (Appendix A) proves valuable 
in controlling this numerical artefact. 

The velfit code can be used to fit simultaneously for 
more than one type of distortion, but we have not ob- 
tained anything us e ful fr om attempting this. However, 
Ivan de Ven fc Fathil (|2009f i suspect both spiral perturba- 
tions and axisymmetric radial inflow in the central parts 
of NGC 1097, and fi t for both using a generaliz ation of the 
method proposed bv lSchoenmakers et al\ (|l997t ). 

It should be noted that velfit is purely a fitting proce- 
dure. It does not, in general, yield a direct estimate of the 
potential responsible for the fitted velocity distortions, ex- 
cept when departures from axial symmetry are mild, as we 
show in the next section. 



3 MILDLY DISTORTED POTENTIALS 

Where the fitted non-circular speeds are a small fraction 
of the circular speed, we can use epicycle theory to relate 
the fitted non-circular velocities to the strength of the non- 
axisymmetric perturbation. We can use the formulae for 
the orbits of test parti cles on near circular orbits in weakly 
barred potentials from lBinnev fc Tremainel (|2008l . hereafter 
BT08)Q Their equation (3.147a, pl90) gives the forced ra- 
dial displacement as a function of time, which can easily 
be differentiated to find the forced radial speed. The radial 
velocity at radius Ro varies sinusoidally with amplitude 



d$ 6 



+ 



2fi$ b 



dR R{Q - fi 6 



m(fio — fit,) 
«q — m 2 (Qo — fib) 2 



(3) 



where <!>(, is the weak non-axisymmetric part of the poten- 
tial that rotates at angular rate fi;,, and fi(i?) & n(R) are 
the usual frequencies of rotation and epicycle motion for 
mildly eccentric orbits. The time der ivative of the equation 
for th e tangential displacement (e.g. ISellwood fc Wilkinson! 
1 19931 . eq. 10b), converted to the same notation and abbrevi- 
ating u — m(fi — fib), gives 



V n 



2fi d$ 6 + 4fij 



dR 



R 



wo 



(4) 



With aspherical halos in mind, we simplify these general 
expressions for the case of a non-rotating potential distor- 
tion in a region where the rotation curve is approximately 
flat. (Other assumptions are possible, if desired.) Since we 
expect the halo to be rotating slowly, we assume fit <C fio, 
we set kq = 2fip for a flat rotation curve, and we assume 
the potential perturbation varies slowly with radius so that 
d&b/dR <C &b/R- With these assumptions, choosing m = 2 
for a bisymmetric distortion and setting i?ofio = V c , equa- 
tions © and @ reduce to 



V 2 ,r 



2$6 

' Vc ' 



Va,, 



3£b 
' Vc ' 



(5) 



Thus if the perturbed velocities are caused by a weak, non- 
rotating, oval distortion to a quasi-logarithmic potential, we 
should expect the perturbed velocity coefficients to be in 
the ratio V^,t — 1.5V^ jT . and, in particular, they should have 
the same sign. Note that we do not expect the perturbed 
coefficients to have this ratio when the oval distortion is 
strong and/or rapidly rotating, such as for bars, or in the 
inner parts where the rotation curve is rising. 

In order to relate to the potential shape, we assume 
a non-axisymmetric potential of the form (c/. eq. 2.71a of 
BT08) 



$(R,e) 



v2 



In 



1 + 



1% 



Ql ■ 2, 

sm l 



(6) 



where Vo sets the velocity scale, R c is the core radius, and 
g<5 is the axis ratio of the potential. In the limit of g$ < 1, 
this is the potential of a midly non-axisymmetric galaxy 
having flat outer rotation curve, of the kind we assumed 
above. Expansion of this potential for small (1 — ?!)/?!>, and 
comparison with the definition of in equations (3.136) 



4 Gas streamlines trace test particle orbits when pressure and 
magnetic forces can be neglected. 
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and (3.143) of BT08, we find $ b = -V?(l - g|_)/(4g|). 
Combining this result with eq. ((5} and equating Vt to V c , 
we finally obtain 

We stress that this formula assumes a mildly distorted, 
slowly rotating potential and a flat rotation curve, and it 
will not yield a reliable estimate in other circumstances. As 
is well known, the density that gives rise to this potential 
is about three times more elongated than the potential, so 
that q p ~ 1 - 3(1 - g$) (BT08, p. 77). 



4 RESULTS 

Here we apply the tool velfit to seve ral galaxies in the 
THINGS sample dWalter et all |200Sl ). The HI Nearby 
Galaxy Survey (THINGS) used the Very Large Array op- 
erated by the National Radio Astronomy ObservatorjQ to 
make spectral observations of the 21cm line emission of neu- 
tral hydrogen in a sample of 34 galaxies. The data are in the 
public domain. We do not reanalyse the entire THINGS 
sample here, but choose a few representative galaxies to il- 
lustrate the advantages of velfit over reswri. 

We selected data with natural weighting, which have 
higher signal-to- noise (S/N) and lower spatial resolution 
than when robust weighting is used, and downloaded maps 
of the intensity-weighted mean velocity. In order to apply a 
S/N cut-off, we also downloaded the data cubes. We deter- 
mined the noise level, a, from parts of channel maps with no 
signal, and discarded velocity measurements from the maps 
for which the peak intensity in any channel < 5a. 

The usual \ 2 function, defined in eq. ([2|, requires an es- 
timate for the error a n in the data at each point in the map. 
Since the neutral hydrogen clouds in a galaxy have typical 
random velocities of 6-12 km s _1 (|Kamphuisl 1 19931 ) . we do 
not adopt the formal (generally much smaller) uncertain- 
ties in the intensity-weighted mean velocities, but instead 
assume a constant uncertainty of a n = 10 km s _1 . We have 
checked that our results are insensitive to values in the range 
<t„ = 10 ±4 km s" 1 . 

We choose the spacing between the rings in our model 
to be the beam width in order that each ring is independent. 
The beam size differs for each galaxy, ranging from 7.41" for 
NGC 2976 to 15.6" for NGC 7793. Since smoothing carries 
with it the danger that it may "smooth away" real features 
in the rotation curve, we do not employ it as a matter of 
routine. However, smoothing is essential to stabilize the fits 
to the data for NGC 2903, where the bar is close to the 
major axis (see J2j) - 

We adopt the distances given in Table 1 o f Walter et all 
(2008) solely for the purpose of marking radii in kpc as well 
as in arcsec in the rotation curve plots. 



5 NRAO is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities Inc. 



4.1 Errors 

Velocity residuals at each pixel have correlations due to fea- 
tures in the data, caused by spiral arms for example, that 
are not part of the fitted model. In addition, a small fraction 
of pixels have large residuals. These properties of the data, 
together with our use of arbitrary uncertainties, imply that 
the values on the \ 2 surface cannot be used to obtain error 
estimates on the parameters of the model. 

We therefore adopt a bootstrap technique to estimate 
the true uncertainties in all fitted quantities, including 
the velocities. This standard non-parametric method (e.g. 
IChernicklll999l ) yields statistical uncertainties in fitted quan- 
tities without making assumptions about the underlying dis- 
tribution of the residuals between the model and the data 
values. The spread of estimated quantities from repeated 
fits to resampled (or pseudo-)data yields an estimate of the 
true uncertainty. We give further details of our procedure in 
Appendix B. 

The error bars on our velocity estimates, and the un- 
certainties in our tabulated values, are the rms variations 
of each quantity from 1000 bootstrap iterations. These sta- 
tistical uncertainties reflect the peculiarities in the distribu- 
tion of the errors in the data values, as well as some of the 
systematic inadequacies of the fitted model, which ignores 
features such as shocks or spiral streaming. They do not in- 
clude possible systematic errors from other sources such as 
beam smearing or pressure support, neither can they reflect 
possible systematic differences between different data sets. 
As we report here, different kinematic data obtained from 
a different instrument, and perhaps arising from a different 
component of the ISM, can yield different estimates of the 
same physical quantities. A disagreement by more than the 
properly estimated errors is an indication either of system- 
atic errors in one or both data sets or the inadequacy of the 
model. 



4.2 NGC 2976: A Case Study 

NGC 2976 is a nearby dwarf galaxy in t he THINGS sampl e 
which was also studied in great detail bv lSimon et al\ (2003). 
Velocity maps have therefore been made using the 21 cm line 
of HI, the optical Ha emission lines, and the 12 CO( J = 1 — > 
0) line. Clear distortions are visible in the velocity maps, 
indicating departures from a simple, coplanar axisymmetric 
flow. 

ISimon et all I^OOSt ) fitted their data with a model that 
combined the usual circular flow pattern with an axisymmet- 
ric radial flow. SS07 found that an oval or bar-like distortion 
could yield an equally good fit to the same data. TBWBK 
find that a tilted ring model fits the HI data, which they 
conclude is consistent with tiny deviations from a round po- 
tential, albeit with a < 30° variation in the inclination and 
a similar change of PA in the inner part of the galaxy. 

Fig. [2ja) shows our fit to the same HI observations used 
by TBWBK when we assume a flat disc plane and fit for 
bisymmetric flows to R — 105" only. The residual velocities 
after subtracting our best fit model are typically less than 5 
km/s, and fewer than 1% of the pixels have residuals larger 
than 20 km/s. The err or bars show ±cr unc ertainties from 
the bootstrap analysis. Ide Blok et all (|2008h report a mild 
warp for R > 125", but we have excluded data at these 
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Figure 2. Our best fit rotation curve and streaming velocities 
derived using velfit on the THINGS data for (a) NGC 2976 and 
(b) NGC 7793. The black line shows the result of an axisymmetric 
fit to all radii and the error bars represent ±<j uncertainties. 



Table 1. Best-fit parameters for NGC 2976 and NGC 7793 
NGC 2976 NGC 7793 



Adopted distance (Mpc) 
Systemic vel. (km s — 1 ) 
Disc inclination i 
Disc PA, <j>' d 
Bar axis (fib 
Projected bar axis (fi' b 



3.6 
1.6 ±0.1 



3.9 

226.9 ±0.1 



61° ± 1° 44° ± 5° 

323° ±1° 292?0±0?4 

-29° ±6° 51° ±8° 

308° ± 3° 334° ± 7° 



large radii because they have low S/N. While our estimate 
of PA (Table[T|) is in Rood agreemen t with that derived from 
the same HI data bv lde Blok et all |200ct ). our value of the 
inclination is somewhat lower than the i = 65° they adopted, 
although their value is estimated from the restricted annular 
range 80" ^ R ^ 110", whereas ours is from a global fit that 
includes bi-symmetric streaming over most of the disc. 

We find (Fig. [2ji) clear non-circular streaming motions 
in the inner galaxy in these HI data that are similar to, 
but of somewhat smaller magnitude than, those found by 
SS07 from the combined Ha and CO data. Even though the 
perturbed velocities are < 10 km s _1 , they are still a large 
enough fraction of the orbit speed that velfit is required; 



notice that the black line, which shows V c estimated in an 
axisymmetric fit such as would be found from reswri, is con- 
sistently lower than Vt. We also find a smaller inclination 
angle than i ~ 64° fitted by SS07. The origin of these differ- 
ences, which are barely consistent within the quoted errors, 
is clearly due to fitting different datasets, with HI data aris- 
ing from a different physical component. 

We estimate the projected orientation of the bar or oval 



to be 



308° ± 3°, which is barely consistent with the 



~ 315° ±4° estimated (from different data) by SS07. The dif- 
ference in angles implies that the bar is slightly farther from 
the major axis, which in turn reduces the difference (Fig.[2K) 
between Vt and the axisymmetric fit, for the reason given in 
section[2] Note, however, that we obtain a projected bar an- 
gle in closer agreement with that estimated by SS07 if we 
fix the galaxy inclination at their estimated value, indicat- 
ing that the non-circular streaming motions in the neutral 
hydrogen arise from forcing by the same non-axisymmetric 
potential. 

The probable reason that the bar streaming velocities 
we find here are weaker than those of SS07 is the signifi- 
cantly less pronounced 'S'-shaped distortions of the velocity 
contours in the HI data (see Fig. 28 of TBWBK) than those 
shown in Fig. 2(a) of SS07. Since the distortions are well- 
resolved, it is unlikely that the differences result from the 
minor difference in spatial resolution of the two datasets: 
the beam-size of the HI data, 7.41", is ~ 50% gr e ater t han 
the 5" smoothing length adopted bv lSimon et all l|2003l ) for 
the combined CO and Ha data used by SS07. (Note that 
ISimon et al\ |2003l ) show that velocity estimates from both 
CO and Ha data, which are almost entirely from HII re- 
gions, are generally in good agreement with each other.) 
We can only speculate as to why the neutral hydrogen has 
milder non-circular motions than either the molecular or ion- 
ized gas; perhaps the HI layer has a greater physical thick- 
ness than the molecular/ionized layer and therefore greater 
"pressure" (in reality larger turbulent speeds), which causes 
a weaker response to a non-axisymmetric potential. Note 
that if the source of the non-axisymmetric potential is a bar 
in the disc of this galaxy, a thick gas layer will feel a weaker 
potential only if the bar is much thinner vertically than is 

the HI layer. 

Since all three models, with radial flows (I Simon et al\ 

l2003h . a twisted disc (TBWBK), or oval streaming (SS07 
and Fig. [2Ji), are adequate fits to the data, there is no sta- 
tistical reason to prefer one over another. However, the oval 
streaming model both avoids the "continuity problem" in- 
herent in radial flow models, and also avoids a strong twist 
in the plane of the inner disc; TBWBK suggest the disc 
plane at R ~ 20" (~ 300 pc) is inclined to the plane of the 
main part of disc at R ~ 1.5 kpc by ~ 30°, while the light 
distribution does not give any indication of such an unusual 
feature. 



4.3 Effects of bar orientation 

The value of Vt in Fig. [2J a), while smaller than obtained 
from different data by SS07, is higher than results from a 
purely axisymmetric flow fit, shown by the black line. The 
bisymmetric fitted Vt is larger in this case because the "bar" 
is oriented such that its principal axis is not far from to the 
major axis of projection for the galaxy, as discussed in sec- 
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tion [2] We therefore also present a case where the bisym- 
metric fit has little effect on Vt. 

Fig. fflb) shows a fit to the THINGS data for 
NGC 7793. As always, w e attempt to fi t only the flat part of 
the disc; in this galaxy, Ide Blok et all (|2008T ) find a mildly 
varying disc inclination over the entire radial range, but the 
PA clearly rises steadily for R > 320", which we tentatively 
interpret as the start of a warp. We fit a bisymmetric flow 
over the inner disk R < 125" only and an axisymmetric 
model to R — 320". The best fit parameters and uncer- 
tainties are given in Table [T] Our best-fit i fc PA are in 
pod agreement with the values estimated by Ide Blok eToH 
2008) from their tilted ring analysis. 

The uncertainty in our estimate of the inclination is, 
however, rather larger than in other cases, possibly because 
of spiral streaming in the outer disc, which is not included 
in our model, or perhaps because the entire disk is warped, 
as suggested by TBWBK. Their explanation seems the more 
likely because the bootstrap values for this parameter have 
a distinctly bimodal distribution symmetrically distributed 
about the best-fit value; conservatively, we estimate the un- 
certainty from the rms spread of the bootstrap iterations. 
The uncertainty in the inclination is reflected in the uncer- 
tainties in the velocities, which therefore seem large relative 
to the smoothly varying means. 

We find clear evidence (Fig.[2p) for non-circular stream- 
ing in the inner parts. Since we find (j>b — 49° , the estimated 
"rotation curve" from the bisymmetric fit is in close agree- 
ment with that from the simple circular flow model, as ex- 
pected from the discussion in section [2] TBWBK describe 
this galaxy as their best candidate for a non-axisymmetric 
potential in the outer parts, but they acknowledge that the 
problem of the apparent changing inclination complicates 
this interpretation. However, the more rapid radial changes 
in i & PA in the inner parts of their fits again appear to 
have masked the non-circular streaming we detect. 
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Figure 3. Our best fit rotation curve and streaming velocities 
derived for NGC 2903 using velfit on (a) the THINGS data and 
(b) the BHaBar data - note the different scales. 



4.4 Strong bars 

TBWBK find non-circular motions that are consistent with 
a round potential in all the THINGS galaxies (their Table 
3), despite the fact their sample contains several galaxies 
that are quite strongly barred. In these cases, however, the 
2fcm line emission from the barred region is generally too 
weak to yield reliable velocity estimates; velocities can be 
measured in limited patches in some cases, e.g. NGC 925, 
but not anywhere in others, e.g. NGC 3627. However, the 
data from NGC 2903 are good enough over almost the en- 
tire barred region (R < 60") to yield velocities above our 
S/N threshold. In this galaxy, TBWBK estimate the me- 
dian value of m = 2 streaming motions in the bar region is 
only ~ 14 km s _1 , and conclude that the mean elongation 
of the overall potential is nevertheless consistent with being 
round. 

The application of velfit to this galaxy presents a sub- 
stantial challenge because the bar is so closely aligned with 
the major axis of projection. As discussed in section [51 the 
velocities Vt, Vi,t & Vi, T become harder to distinguish as 
\4>b | — > 0, and velfit can return unphysically large values for 
all (S> 1000 km s _1 in magnitude), even though they still 
fit the data well when combined as eq. (JTJ. We overcome 
this problem to a large extent by applying a very small 



Table 2. Best-fit parameters for NGC 2903 

THINGS data BHaBar data 



Adopted distance (Mpc) 
Systemic vel. (km s -1 ) 
Disc inclination i 
Disc PA, <j>' d 
Bar axis (pf, 
Projected bar axis (j>' b 
Smoothing penalty, A 



8.9 

549.9 ± 0.3 
64° ± 1° 
201?5±0?5 
6° ± 14° 
204° ± 6° 
1.7 x 10~ 4 



554.1 ±0.5 
66° ± 3° 

204° ± 1° 
-12° ±8° 

199° ±4° 
6.6 x 10" 5 



smoothing penalty, as described in Appendix A. Fig[3Ja) 
shows the results from velfit applied to the THINGS data 
for R < 400"; the parameters of this fit are given in Table [2] 
for A = 1.7 x 10~ 4 , but results are insensitive to variations 
of a factor of a few about this value 

The smoothing penalty successfully eliminates absurd 
velocities in the bar region obtained from the bootstrap anal- 
ysis, but at the cost of introducing a strong bias against find- 
ing \cj>b\ < 2°. The reason appears to be a ridge in the \ 2 sur- 
face as 0b — > caused by the smoothing penalty, which dis- 
favours the wildly varying velocities that would achieve the 



6 This apparently bizarre choice is A 
200 km s — 1 ; see Appendix A. 



0.002/.4 for Vtyp 
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smallest residuals. We therefore prefer a very small smooth- 
ing penalty since larger values widen the range of disfavoured 
bar angles even though they further reduce the scatter of 
fitted velocities from the bootstrap iterations. On the other 
hand, increasing the smoothing penalty 50-fold, leads to an 
almost linear rise in Vt, even for the bisymmetric fit, which 
illustrates the perils of oversmoothing. 

Our best-fit estimate of the bar angle is 4>b ~ 6°. It 
should be noted that neither this angle, nor the best fit val- 
ues of Vt change significantly when we eliminate the smooth- 
ing penalty altogether. 

Our fitted valu es of i fc PA a r e in g ood agreement with 
those estimated bv lde Blok et al] (|2008l ). and their rotation 
velocities are in good agreement with our axisymmetric fit, 
the black line in Fig.[3ja), that rises in a quasi-linear fashion 
from the origin. Naturally, they find wildly varying values of 
both i and PA in their innermost few rings, whereas velfit re- 
quires a flat plane and consequently our bisymmetric fit (red 
line) finds different velocities in this region. We find some- 
what larger non-circular speeds in the bar region than those 
estimated by TBWBK, although the uncertainties in our es- 
timates also are large enough that a round potential could 
not be excluded. Thus, the surprisingly weak bar streaming 
motions are not merely an artefact of reswri. 

Since the beam width of the HI data we fit is ~ 15" 
(equal to our ring spacing in velfit) and the bar semi-axis is 
about four beam widths, it is likely that the HI data is not 
fully able to resolve the bar flowQ Fortunately, the velocity 
field of this galaxy has also been mappe d at higher spatia l 
resolution in the 12 CO(J = 1 -> 0) line dHelfer et fflZj|2003h 
and i n Ha using a Fabry-Perot instrument ( Hernandez et all 
120051 ). We here show fits to the Ha data (kindly made avail- 
able by Olivier Hernandez) since they extend to larger radii 
than do the CO data, albeit with lower spatial resolution 
(4.8"). 

Fig. 0Ja) shows a 2MASS0 H-band image of the galaxy, 
together with (b) the velocity map from the BHaBar sur- 
vey, (c) our best fit model and (d) & (e) residuals from an 
axisymmetric and full bar flow fits. 

Table [2] and Fig.[3jb) give the results we obtained from 
by applying velfit to the Ha data for NGC 2903 - note that 
the radial scales in panels (a) & (b) of Fig. ^differ. We used 
a similar smoothing (A = 6 x 10 -5 in this case) in order to 
eliminate absurdly large velocities in the bootstrap analysis, 
which introduces a bias, as before, against bar angles close 
to zero; since this bias causes a strong skewness in the dis- 
tribution of bar angles, we estimated the uncertainty in this 
quantity from only those values more negative than the best 
fit value. The parameters of the fit to these data are gen- 
erally in good agreement with those from the HI data, and 
the projected bar angles differ within their uncertainties. 

However, the Ha data show a much more pronounced 



7 Robust weighting of HI data yields velocity maps with higher 
spatial resolution but lower S/N. For NGC 2903, however, these 
data from the bar region have too low S/N to allow meaningful 
fits. 

8 Atlas Image obtained as part of the Two Micron All Sky Survey 
(2MASS), a joint project of the University of Massachusetts and 
the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 



non-circular flow pattern within the bar region (R < 60"), 
with perturbed velocities almost as large as Vt, which in 
turn significantly exceeds the estimated circular speed from 
an axisymmetric fit to the same data (black line) and the 
derived Vt from the HI data. The fitted velocities from the 
two datasets are in reasonable agreement outside the bar re- 
gion, for as far as the optical data extend. The uncertainties 
in the bar region are still large, and significantly larger than 
the point to point variation in the best fit suggesting that 
slightly more aggressive smoothing could be warranted. 

Thus it is clear that inadequacies in the HI data are the 
reason TBWBK concluded that non-circular motions within 
the bar were small. While the beam size of the THINGS 
data is not fully adequate to resolve the bar flow, the general 
paucity of neutral gas in the barred regions of other galaxies 
suggests it is also likely that HI is simply not a faithful tracer 
of the bar flow. Data from a different component of the ISM 
having better spatial resolution and fuller spatial coverage 
do reveal a pronounced non-circular streaming pattern, as 
expected for this strongly barred galaxy. 



4.5 Searching for Aspherical Halos 

We here attempt to constrain the shapes of dark matter 
halos by searching for non-circular streaming motions in the 
outer discs of the THINGS sample. Since the current version 
of velfit assumes the galaxy plane to be flat, it cannot be used 
in warped regions, which generally arise outside the visible 
disc. We are therefore restricted to just two galaxies in the 
sample, NGC 3198 & NGC 2403, fo r which the ex t ended HI 
disc is known from the analysis of Ide Blok et al\ |2008T l to 
be approximately coplanar with the inner disc. 

Even though these two galaxies are not strongly warped, 
the analysis of TBWBK, which allows changes in PA and i 
from ring to ring, may underestimate the ellipticity of the 
dark matter halos in the disc plane. However, the main ad- 
vantage of velfit over reswri in these circumstances is that it 
searches for a bisymmetric distortion that is coherent over 
a wide range of radii and could, in principle, detect very 
mild distortions that might be masked by various sources of 
noise, such as turbulence and local spiral streaming. Since 
it is a more sensitive probe of halo shapes, it should either 
detect mild distortions, if they are present, or place a tighter 
lower bound on the axis ratio of the potential. 



4.5.1 NGC 3198 

Fig. [5] shows the residual map to R = 456" when a flat, 
axisymmetric model is subtracted from the THINGS data 
for NGC 3198. The residual velocities are generally small, 
peaking at ±20 km s~\ which is consistent with the small 
variat ions in i & PA for r > 200" reported bv lde Blok et al] 
(2008). However, the residual pattern reveals clear indica- 
tions of mild spiral arm streaming, ev en far outside the op- 
tical disc (-R25 — 255" in the B-band. Ide Vaucouleurs et all 

In order to search for a possible mildly non- 
axisymmetric halo, we tried fitting a bisymmetric model, 
with no spirality, to the outer disc. Such a model may be 
able to identify a weak bar flow that could be buried in the 
spiral noise. The parameters of our best fit model, which 
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Figure 4. (a) An H-b and image of NGC 2903 from 2MASS llSkrutskie et aLll2006h . (b) The velocity map made using the Ha line 
l|Hernandez et alll2005h . (c) Our best fit model with a bar flow, (d) Residuals after subtracting the best fit axisymmetric model - note 
the large values in the bar region, (e) Residuals after subtracting our best fit bi-symmetric flow model. The velocity scales to the right 
are in km s — 1 ; all images have the same orientation and spatial scale, indicated by the V scale bar. 



Table 3. Best-fit parameters for NGC 3198 and NGC 2403 
NGC 3198 NGC 2403 



Adopted distance (Mpc) 
Systemic vel. (km s _1 ) 
Disc inclination i 
Disc PA, <j>' d 
Bar axis <fib 

Projected "bar" axis <j>' b 



13.8 
660.2 ±0.3 
70?4±0?3 
216? 2 ± 0?5 
46° ± 14° 
236° ± 5° 



3.2 

133.5 ±0.3 
64° ± 1° 
123?8 ± 0?6 

296° ± 12° 
82° ± 5° 



includes a bisymmetric flow for R > 220", are listed in Ta- 
ble [3] Our estimated values of V S v S , i, fc PA are in e xcellent 
agreement with those given bv lde Blok et al. | (120081 1. 

Fig. [6] shows that the best fit non-axisymmetric veloci- 
ties are no larger than ~ 6% of the circular speed and vary 
slowly with radius. If the perturbed velocities are caused 
by a slowly-rotating, mild oval distortion of the halo in the 
outer parts where the rotation curve is approximately flat, 
we should observe V-z.t — \SV%, r (eq. [S]). In fact, the coeffi- 
cients are not in this predicted ratio at any radius and even 
the signs differ over the inner part of the fitted range sug- 
gesting a different origin for the perturbed velocities, such as 
spiral arm streaming for which a fixed axis and slow rotation 
for the perturbation are inappropriate assumptions. Our fit- 
ted values of 06, V2,t(R), and V2, r (R) are merely those that 
achieve the largest reduction in \ 2 from the residual pattern 
shown in Fig. [5] In particular, the PA of the fitted "bar" 
varies with radius as we sub-divide the fitted region, as one 
would expect if the fit is picking up different fragments of 
the spirals. 

Thus we are unable to identify a velocity pattern in- 
dicative of a non-axisymmetric halo. Whatever possible halo 
distortion may be present, it is clearly still weaker than the 
mild spiral features we do detect. We therefore concur with 



TBWBK that the data from NGC 3198 are consistent with 
this galaxy living in a perfectly round halo. 

Even though we do not have a firm detection of a bar- 
like distortion in the halo, we can use the results in Fig.[6]to 
place a lower bound on its ellipticity by making a rough esti- 
mate of the maximum halo distortion that could be masked 
by the spiral contribution. We therefore suppose that the 
disturbed velocities are due to a combination of a slightly 
elongated halo and a spiral disturbance in the disc plane. 
In this situation, the variation of the total amplitude of the 
velocity distortions will be modulated by the changing angle 
between the two as the phase of the spiral changes with ra- 
dius; it will be a maximum when the distortions are aligned 
and a minimum when they are perpendicular. (This is a con- 
servative assumption, since amplitude and phase variations 
could also be caused by the changing distance from a reso- 
nance in just a single spiral, with no second component of 
the disturbance potential.) 

For NGC 3198, the combination A to t = (Vl,, + V 2 2 t ) 1/2 
varies from a minimum of A to t = 1.25 km s _1 at 296" to 
-Atot = 11-10 km s" 1 at 444". If this spans a full range of 
possible phase differences (0 to 7r/2), and the spiral and halo 
contributions are separately roughly constant, then -Atot = 
A s + Ah at maximum, and -Atot = A s — Ah at minimum. 
Solving for the separate spiral and halo amplitudes, we find 
A s ~ 6.3 km s _1 and Ah — 5.0 km s _1 . Since V% t t = 1.5V2, r 
for a mild, non-rotating distortion with a flat rotation curve 
(eq. [5]), we expect Vi, T — Ah/1.8. Using the values V2,r — 
5/1.8 km s" 1 and Vt = 140 km s _1 in formula ©, we find 
g$ > 0.98 and q p > 0.94 as our estimated lower limits on 
the axis ratios of the potential and density of the halo in 
NGC 3198. 
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Figure 5. Residuals after fitting an axisymmetric model to the 
THINGS data for NGC 3198. The outer ellipse has a semi-major 
axis of r = 456" and velocities in the colour bar are in km s _1 . 
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Figure 6. Results from velfit using THINGS data for NGC 3198. 
The black line shows the result of an axisymmetric fit to all radii, 
while the coloured points show Vt(R) (red), V2,t(R) (green), &: 
V2,r(R) (blue) from a bisymmetric fit restricted to the range 
240" ^ R ^ 460" with an axisymmetric model fitted to other 
radii. 



4.5.2 NGC 2403 

Fig. [7] shows that a simple, flat, axisymmetric model is a 
good fit to the THINGS data for NGC 2403, c onsistent with 
the t iny variations in i & PA reported by (|de Blok et all 
2008). The small residuals are somewhat less indicative of 




Figure 7. Residuals after fitting an axisymmetric model to the 
THINGS data for NGC 2403. The outer ellipse has a semi-major 
axis of r = 844" and velocities in the colour bar are in km s — 1 . 
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Figure 8. Same as for Fig. [6] but for NGC 2403. 



spiral streaming than in NGC 3198, and show hints of a more 
global N-S anti-sy mmetry at larger projected radii (R25 — 
656" in the B-band lde Vaucouleurs et aZ.lll99lh . 

Our fitted values of V ByB , i & PA from a model that in- 
cludes bisymmetric streaming velocities for R > 450" , listed 
in Table [3l are in ex cellent agreement with those found by 
Ide Blok eioH (|200ot ). Fig. [8]shows that the perturbed veloc- 
ities are < 10 km/s, but the radial component is generally 
the larger, which is inconsistent with a non-rotating, bar- 
like distortion. Thus we again concur with TBWBK that 
the THINGS data on this galaxy are consistent with it liv- 
ing in a perfectly round halo. 

While the perturbed velocities are not of the form ex- 
pected for a mild, non-rotating halo distortion, we pro- 
ceed as for NGC 3198. For NGC 2403, we find Act varies 
from a maximum of Aot = 11.86 km s^ 1 at 510" to 
Atot = 2.12 km s _1 at 809". Again assuming this spans 
a full range of possible phase differences and the spiral 
and halo contributions are separately roughly constant, then 
A a ~ 7.0 km s _1 and Ah — 4.9 km s _1 . Using the values 
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V 2 , r ^ 4.9/1.8 km s" 1 and Vt = 120 km s" 1 in formula ©, 
we again find q<s> > 0.98 and q p > 0.94 as our estimated 
lower limits on the axis ratios of the potential and density 
of the halo in NGC 2403. 



5 CONCLUSIONS 

ISpekkens fc Sellwoodl (2007) devised the software tool velfit, 
for fitting non-circular streaming flows in galaxies caused 
by non-spiral like distortions. We have shown here that 
it is altogether superior to the commonly used reswri 
ijSchoenmakers et aZ.I Il997l ). This is because it provides a 
correction for the systematic error in the azimuthally aver- 
aged orbital speed, it can fit strong distortions, it is more 
sensitive to radially coherent disturbances, it readily allows 
radial smoothing of the fitted velocities, and the estimated 
disturbed velocities are more easily related to the potential 
distortion. We also argue that while the current version of 
velfit assumes a flat plane for the inclined disc, this is no 
more restrictive, since tilts of individual rings in the bright 
inner part of the disk, as reswri allows, are dangerous and 
can lead to severe underestimation of the distorted veloci- 
ties. This paper illustrates these advantages for a number of 
galaxies. 

We show that the THINGS data (jWalter et al.ll2008T l 
for NGC 2976 can indeed be fitted by an inner bar-like dis- 
tortion, albeit somewhat milder than that found by SS07 
from Ha+CO data. We argue that a ba r flow is more natu - 
ral than either the radial flo w fitted by Simon et al\ (I2003T I 
or the twisted disk model of iTrachternach et al\ (|2008T l. As 
shown in SS07, the improved mean orbital speed of the gas 
estimated from velfit is a fairer estimate than the simple 
mean fitted by tools such as reswri. The difference, which 
arises from the bias to the velocity on the major axis, can 
be of either sign depending on the orientation of the bar to 
the major axis of projection. 

We also show that neutral hydrogen observations are 
not well suited to tracing gas dynamics in strongly barred 
potentials. The neutral hydrogen generally has a low col- 
umn density in the barred region, and the velocity maps are 
generally too noisy or sparsely sampled to yield a clear indi- 
cation of bar flow. Smoothing to lower spatial resolution 
improves signal-to-noise, and yields an almost fully sam- 
pled velocity map throughout the bar region of NGC 2903. 
We identify an oval flow pattern in these data of about the 
right physical size and with signifiant streaming velocities, 
but the large uncertainties imply they are also consistent 
with a round potential. However, a strong bar flow is un- 
a mbiguously detect e d in our fit to the Ha velocity map 
of lHernandez et al\ (|2005l ). While the superior spatial res- 
olution of the optical data is clearly important, the gen- 
erally patchy and faint emission from neutral hydrogen in 
the barred regions of this and other barred galaxies in the 
THINGS sample suggest that neutral hydrogen is simply a 
poor tracer of bar flows. 

Our analysis of the THINGS data for NGC 3198 & 
NGC 2403 reaches a similar conclusion to that of TBW BK: 
that t he h alos of these t wo ga laxies are close to round. |jod 
(2000) and lBailin etaH ((2007) show that the self-consistent 
response of the disc can mask a large part of the distortion 
in the halo, but only when the disc is massive. In these 



two cases, the outer gas disc probably has little mass and 
therefore could not hide a more substantial halo distortion. 

However, it is hard to confront this result with the pre- 
dictions of LCDM (cited in the introduction): not only is 
it based on just two galaxies, but it is possible these two 
galaxies are unrepresentative perhaps because only galax- 
ies with unusually round halos could host an extensive disc 
of gas that is not warped! Clearly, measurements of halo 
shapes in a representative galaxy sample will require a tool 
that can reliably measure potential distortions in warped 
discs; we leave development of suc h a capability for future 
work. Other statistical ap proaches (|Franx fc de Zeeuwlll992l ; 
ITrachternach et qZ.I l2009h do, however, place some reason- 
ably tight constraints on halo shapes. 

Thus, while we confirm the conclusion of TBWBK that 
many galaxies in the THINGS sample have at most minor 
departures from axial symmetry, velfit reveals that mild bars 
are present in NGC 2976 and NGC 7793. We also show that 
a pronounced non-axisymmetric flow is revealed more clearly 
in other data for the strongly barred galaxy NGC 2903. 

Another valuable application for velfit will be a prelim- 
inary analysis of the velocity maps of strongly barred galax- 
ies. It would be very helpful to obtain a clear indication of 
whether the flow pattern is, or is not, well enough sampled 
and sufficiently regular to yield a result, before embarking 
on laborious mass modeling by the method described by 
IWeiner et al. (2001). Furthermore, such a study needs ap- 
proximate axisymmetric mass models, and the estimates of 
Vt from velfit will be more useful than the "circular speed" 
estimated from other less powerful tools. 
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APPENDIX A: APPLYING A SMOOTHING 
PENALTY 

The usual x 2 function is defined in eq. but we can min- 
imize a new function if we wish 



x ' 2 = x 2 + A J2 A t^- 1 - m + 14 + 1 ] 2 



(Al) 



which adds a penalty for large second differences between 
the tabulated values Vk, which are assumed to be equally 
spaced. The constant A, defined below, has dimensions of 
inverse square velocity in order to fulfil the requirement that 
the smoothing penalty be dimensionless. The value of the 
smoothing parameter, A, which can be set independently 
for the axisymmetric and non-axisymmetric terms, can be 
adjusted as desired; small values have slight effect, while 
large values make the profile very smooth. Over-smoothing 
not only increases the value of \ 2 > but may reduce the inner 
slope and smooth other real features in the rotation curve. 

Since the magnitude of the second velocity difference (in 
square brackets) varies as the inverse square of the ring spac- 
ing, AR = Rm^/N R , we choose A = (iCax/Ai^/Kyp = 
Nft/Viyp, which ensures that the smoothing constraint is 
unaffected when the ring spacing is changed. Here, Vt yp is a 
constant that is a rough estimate of a typical orbital speed 
in the disc. 

The smoothness constraint is applied only in the matrix 
that is solved to find the optimal values for Vk (see SS07), 
and adds extra terms as follows: 



OX 2 
dV 



dVi 



(A2) 



+2A A (Vj-2 - AVj-i + 6Vj - W j+ i + V J+2 ) . 



APPENDIX B: BOOTSTRAP ESTIMATION OF 
ERRORS 

The bootstrap uses the spread of estimated quantities from 
repeated fits to resampled (or pseudo-)data to yield a non- 
parametric estimate of the true uncertainty in each quantity. 
To construct one realization of pseudo-data, we add to the 
predicted velocity from the best-fit model at every pixel a 
residual velocity from some other pixel chosen at random 
from residuals between the best fit model and the real data. 
Completely random resampling of residuals assumes that the 
residuals are uncorrelated, whereas inspection of our residual 
images generally reveals coherent patterns of residual veloci- 
ties due to features, such as other forms of non-axisymmetric 
streaming and large-scale turbulence, which are not included 
in our fitted model. Pseudo-data constructed by fully ran- 
dom interchange eliminates correlations between the residu- 
als and generally leads to unrealistically small uncertainties. 
Thus we require a scheme that will reproduce appropriately 
correlated errors at random at each iteration of the boot- 
strap^ 

ISpekkens fc Sellwood! l|2007h adopted constant residu- 
als over small patches of the image, which worked well for 
the small galaxy NGC 2976 where most of the correlations 
appeared to arise from turbulence. But that scheme cannot 
capture the patterns of residuals that arise from spirals and 
other non-axisymmetric distortions that are clear features 
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in the residuals for larger galaxies. We therefore adopt a 
different approach here. 

Since the dominant residual correlations appear over 
azimuthally extended regions in the disc plane, we manipu- 
late the actual residual pattern as follows: we deproject it to 
face on and then, at every iteration, we rotate the residuals 
through a random angle. We also shift the residual pattern 
outwards by adding a constant to the radius of every pixel 
and subtract the maximum radius from those pixels that are 
shifted outside the map so that they fill the hole created in 
the centre by the outward shift. The constant used to shift 
in radius is a randomly-chosen fraction of the radius of the 
map. We then reproject the new residual pattern, and as- 
sign a residual velocity at each pixel in the pseudo-data from 
that of the nearest pixel in the scrambled residual image. 

When we fit for an axisymmetric model in part of the 
galaxy and include a non-axisymmetric perturbation over a 
limited radial region, we scramble the residuals within each 
of these parts separately, and do not interchange residuals 
between the separate regions, since the appearance of resid- 
ual patterns in the two regions can differ quite markedly. 
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